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Abstract 



We analytically work out the long-term rates of change of the six osculating Kcplcrian 
orbital elements of a test particle acted upon by the Lorentz-violating gravitomagnetic ac- 
celeration due to a static body, as predicted by the Standard Model Extension (SME). 
Wc neither restrict to any specific spatial orientation for the symmetry-violating vector 
s = {— s"^, — s°^, — s*'^} nor make a priori simplifying assumptions concerning the orbital 
configuration of the perturbed test particle. Thus, our results are quite general, and can be 
applied for sensitivity analyses to a variety of specific astronomical and astrophysical scenar- 
1/^ ' ios. We find that, apart from the semimajor axis a, all the other orbital elements undergo 

OO , non- vanishing secular variations. By comparing our results to the latest determinations of 

the supplementary advances of the perihelia of some planets of the solar system wc prclim- 
Cn ! inarily obtain s^ = (0.9 ± 1.5) x IQ-^, Sy = (-4 ± 6) x 10-^ s^ = (0.3 ± 1) x IQ-^. Bounds 

O ■ from the terrestrial LAGEOS and LAGEOS II satellites are of the order of s - 10'^ - 10"''. 

(N 

PACS: 04.80.-y; 04.80.Cc; 04.50.Kd; ll.30.Cp; 95.10.Km 
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There are a variety of theoretical schemes predicting violations of Lorentz symmetry. Attempts 
to quantize gravity have resulted in theories that allow for it [I] . Certain string theories envisage 
the possibility of a spontaneous breaking of it [2]. The possibility of detecting experimental 
signatures of breakings of Lorentz symmetry [3l [1] from an underlying unified theory at the 
Planck scale has recently raised interest in Lorentz violation; for a recent review of modern 
tests of Lorentz invariance see, e.g., [Ij. 

The Standard Model Extension (SME) is a theoretical framework which allows for generic 
violations of the Lorentz symmetry for both gravity and electromagnetism O O [71 [8]. In 
general, there are 20 coefficients for Lorentz violation in the gravitational sector; by assuming 
spontaneous Lorentz-symmetry breaking, the main effects in the weak-field approximation are 
accounted for by the traceless coefficients s'^'^ [9] containing 9 independent quantities. 

According to Bailey [10], in the weak field and slow motion approximation, a test particle 
moving with velocity v at distance r from a central, static body of mass M experiences a 



Lorentz-violating gravitomagnetic acceleration 



Agm = (^) X Bg, (1) 

where 

^G = -^:^ (^ X ^) ' (2) 

with 

s = -s°^ j = 1,2,3; (3) 

G is the Newtonian constant of gravitation, while c is the speed of light in vacuum. Constraints 
from Lunar Laser Ranging (LLR) jllj are |12j 

s"^ = (-0.8±Ll) X 10-^ (4) 

s°2 = (-5.2 ± 4.8) X 10"^, (5) 

so that 

Smax < 1 - 2 X 10-^ (6) 

A general account of the present-day bounds on all the SME Lorentz violating coefficients can 
be found in [T3] . 

From simple dimensional considerations, planetary orbital precessions in the field of the 
Sun, if present, would be as large as 



^ 



Ag 2GMs 1.77 X 10^ m2 s"^ 



< — - = 12 ' (7) 



V cv^ 



where ^ denotes a generic orbital element, and we used eq. ([6]). The present-day level of 
accuracy in constraining the secular rates of orbital changes of the planets of the solar system 
are listed in Table [1] and Table [21 In Table [3] we show tentative upper bounds on putative 
precessional effects according to eq. ([6]) for all the planets of the solar system. A tension between 
the LLR bounds [12] of eq. Q-eq. ([5]) and the resulting expected planetary precessions of Table 
[3] exist, especially for the inner planets. 

Thus, a more detailed analysis is worthwhile. Actually, a mere order-of-magnitude analysis 
based just on eq. ([7]) and Table[3]would be insufficient to draw meaningful conclusions. Lideed, 
exact calculations of the secular variations of all the osculating Keplerian orbital elements caused 
by eq. ([T]) must be performed, with, e.g., standard perturbative techniques, in order to check if it 
really induces averaged non-zero orbital changes. Moreover, also in such potentially favorable 
case caution is still in order. Indeed, it may well happen, in principle, that the resulting 
analytical expressions retain multiplicative factors l/e ,k = 1,2,3,... or e ,k = 1,2,3,... 
which would notably alter the size of the found non-zero secular rates with respect to the 
expected values according to eq. ([7]). Thus, in Section [21 we will analytically work out the 
long-term rates of change of all the osculating Keplerian orbital elements of a test particle 
acted upon by eq. ([T]) by using the Gauss perturbative scheme [17]. 

2 Analytical calculation of the orbital precessions 

In the following, we will not make any particular assumption about the orientation of s in 
space, so that we will adopt a generic reference frame centered in M, with respect to which 



Table 1: Uncertainties in the rates of change of the semimajor axis a, the eccentricity e, the 
inchnation / to the mean echptic at J2000.0, the longitude of the ascending node H, the 
longitude of perihelion w = Q + to, the mean motion n^ = VGMa~^ and the mean longitude 
A = n + uj + A4 of the planets of the solar system; to is the argument of perihelion, while M 
is the mean anomaly. They were inferred by us by taking the ratios of the formal, statistical 
errors in Table 3 of [14j, all rescaled by a factor 10, to the data time span AT = 93 yr (1913- 
2006) of the EPM2006 ephemerides used by Pitjeva [H]. The figures for rib account also for the 
uncertainty acM = 10 km^ s~^ in the Sun's gravitational parameter GM retrieved from |15j . 
The results for Saturn are relatively inaccurate with respect to those of the inner planets since 
radiotechnical data from Cassini were not yet processed when Table 3 of [H] was produced. 
Here mas cty~^ stands for milliarcseconds per century. 
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Mercury 


3.6 




4 X 10"^ 


14.8 


121 


5.4 


50.6 


3.7 


Venus 


2.3 




2 X 10~i° 


0.3 


9.9 


5.7 


10.2 


0.3 


Earth 


1.5 




5 X 10-11 


— 


— 


0.6 


5.2 


— 


Mars 


2.8 




5 X 10-11 


0.03 


0.8 


0.1 


2.8 


0.02 


Jupiter 


6612.9 




2 X 10-*^ 


23.4 


1138.9 


79.8 


131.2 


15.7 


Saturn 


45763.4 




1 X lO-'^ 


42.7 


806.7 


778.2 


196.9 


37.6 


Uranus 


433269.0 




3 X 10-^ 


66.2 


3671.9 


836.5 


323.7 


80.7 


Neptune 


4.9818 X 


10« 


8 X 10-^ 


75.3 


2284.0 


20249.1 


1210.7 


263.2 


Pluto 


3.66961 > 


:107 


3 X 10-6 


167.4 


305.4 


2713.4 


4501.6 


468.3 



Table 2: Supplementary advances of perihelia and nodes of some planets of the solar system 
estimated by [16] with the INPOPlOa ephemerides. Data from Messenger and Cassini were 
included for Mercury and Saturn. The reference {x, y} plane is the mean Earth's equator at 
J2000.0. 



^(^ ^^ 



Mercury 


1.4 ±1.8 


0.4 ±0.6 


Venus 


0.2 ±1.5 


0.2 ±1.5 


Earth 


0.0 ±0.9 


-0.2 ±0.9 


Mars 


-0.05 ±0.13 


-0.04 ±0.15 


Jupiter 


-40 ± 42 


-41 ± 42 


Saturn 


-0.1 ±0.4 


0.15 ±0.65 



^ = {sx,Sy,Sz}- Moreover, we will not make any definite choice about the reference {x,y} 
plane, so that our results can be applied to a variety of specific astronomical and astrophysical 
scenarios. 

Thqj radial, transverse and normal components Aji,AT,A]\f of eq. ([I]), evaluated onto the 



^The radial unit vector R is directed from A/ to the test particle, the transverse unit vector t lies in the 
osculating orbital plane and is perpendicular to R, while the normal unit vector N is perpendicular to the orbital 
plane and is directed along the osculating orbital angular momentum L. Such three unit vectors constitute a 



Table 3: Orders of magnitude of the maximum values 



precessions due to eq. ([T|), with s = Smax 
[T] and Table [2l 



^ 



of putative planetary orbital 



2 X 10 ^. Cfr. with the empirical bounds in Table 
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mas 
cty 



^ 



U- 

Vcty 



1.7 X 10"*^ 
4.7 X 10""^ 

2.5 X 10""^ 
1.1 X lO"'^ 
9.3 X 10-^ 
2.7 X 10~9 
6.7 X 10~^° 
2.7 X 10~^° 

1.6 X 10-1° 



Mercury 


343.70 


Venus 


98.43 


Earth 


51.50 


Mars 


22.18 


Jupiter 


1.92 


Saturn 


0.56 


Uranus 


0.14 


Neptune 


0.05 


Pluto 


0.03 



unperturbed Keplerian ellipse, are 

-2GMnb(l + ecos/)\ 
^ = 7, ^^2 ^^' 



(8) 



•^Tl = Sz cos u sin / + cos / cos u {sy cos fi — Sx sin Q) — sin u [sx cos /? + Sy sin Q) 



(9) 



At 



2eGMnb(l + ecos/) sin/ 



ac(l — e^) 



2\5/2 



^r, 



(10) 



^7- = Sz COS w sin / + cos / cos u {sy cos H — Sx sin H) — sin it (s^ cos H + Sy sin J?) , (11) 



. 2eGMnb (1 + e cos /)^ sin / . 

^N = z ^^57^ -^^' 

ac(l — e"^) ' 



(12) 



•^Af = Sz COS / + sin / {sx sin Q — Sy cos J?) , (13) 

where [IT] u = / + cj is the argument of latitude defined in terms of the true anomaly / and 
the argument of pericenter a;. Notice that At = A^ = for circular orbits, i.e. for e — )• 0. 
Inserting eq. ([Sj-eq. ()12p in the right-hand-sides of the Gauss equations for the variation of 
the elements and averaging them over one orbital revolution yield the following non-vanishing 
rates of changes of the Keplerian orbital elements 



da 
IE 



0, 



(14) 



right-handed orthonomial basis comoving with the test particle. 



de\ _ 2GMV1 - e2 



at I ^^2 



1 + Vl - e^ 



(15) 



£^ = cosw (sa; COS i7 + Sy sin /?) + sinc;; [s^sin/ + cos/ (sy cos Q — s^ sin J?)] , (16) 



dl 
'dt 



2GM (l - VT^^ 



■:?:, 



X = sin w [s^ cos / + sin / [sx sin Q — Sy cos i7)] , 



(17) 
(18) 



dfl 
'dt 



2GM (l-VT^^' 



ca^eVl — e^ 



■AA, 



AA = cos a; esc / [sz cos / + sin / {sx sin Q — Sy cos i?)] , 



droX _ 2GM 

~di/ ~~ca2e3(l-e2) 



^, 



(19) 
(20) 
(21) 



"P = (l — e^) ( 1 — V 1 — e^ j sinc^j (s^; cos fi + Sy sm. Q) + 



+ cosuj{sy (? f-l + e^ + \/l -e^") + 

+ (-1 + ^1-62-6^ (-2 + e^ + 2V1 - 62)) cos / 



cos i7+ 



+ (l — e ) (1 — vl~e2js2 sin/ — e v 1 — e2sxsin i? — v 1 — 6253. cos /sin J7+ 



+ 2e^ \/l - e2sa; cos / sin J7 + (-1 + e^) Sx [-e^ + (-1 + e^) cos/] sin/2- 



e v 1 — e2s^cos 



/ tan (~)+e (l — e)s^ cos / tan I — 



(22) 



dM\ _ 1GM (1 - e2) 



dt 



cd^e ( 1 — e2 + Vl — &' 



X, 



(23) 



11 = Sz COS w sin / + cos / cos a; (sj^ cos fl — Sx sin /?) — sin a; (s^- cos 12 + Sj^ sin J?) . (24) 



The results of eq. (|14p -eq. (j23p are exact in the sense that no a priori simpUfying assumptions 
on either e and / have been assumed. By expanding in powers of e, it turns out that, in the 
hmit of small eccentricities, the inclination and node precessions are of order 0{e), while the 
precessions of the perihelion and the mean anomaly are not defined when e — )• since they are 
of order 0{e~^); on the other hand, the precession of the mean longitude X = lu + J^ is of 
order 0{e). For other computation of the orbital effects, see [9]. The generality of eq. (fH|) -eq. 
(|23p allows one to use them for sensitivity analyses in a variety of specific astronomical and 
astrophysical scenarios, where different {x, y} reference planes and orbital configurations of test 
particles appear. 

3 Results and conclusions 

By using eq. (|2ip for Mercury, Venus and the Earth and the figures of Table [2] for their 
supplementary advances of perihelia, we solve for Sx,Sy, Sz and obtain 

Sx = -s°^ = (0.9 ± 1.5) X 10~^ (25) 

Sy = -s°^ = (-4 ± 6) X 10"^ (26) 

Sz = -s°^ = (0.3 ± 1) X 10"^ (27) 

The forthcoming analysis of more data from the ongoing MErcury Surface, Space ENvironment, 
GEochemistry, and Ranging (MESSENGER) mission [18] to Mercury should allow one to im- 
prove such bounds in a near future. In view of that fact that, according to eq. (|14p -eq. (|23p . 
the strongest signals occur for the closest particles to the central body, it would be desirable 
that supplementary advances of all the Keplerian orbital elements of Mercury will be produced 
in future ephemerides. 

We stress that eq. ()25p -eq. ()27p should be regarded as a sort of predicted parameter 
sensitivity since, actually, Fienga et al. [16] did not model any Lorentz-violating terms in 
the INPOPOlOa ephemerides. Actually, in more refined analyses, it would be required, in 
principle, to explicitly include eq. ([T|) in the dynamical force models usually fitted to planetary 
observations, and solve for the s'^-', j = 1,2,3 parameters in a dedicated global solution obtained 
by reprocessing the entire data set with such modified softwares. The use of a similar approach, 
including also LLR data, was envisaged in [12]; see also the considerations by Nordtvedt [19] 
in a different context. 

The Moon yields us a benchmark for testing the degree of reliability of our approach. Indeed, 
from the small eccentricity limit of eq. (|2ip . it can naively be posed 

GMs , , 

w —, 28 

from which one can infer 

'<^,i^. (29) 

Since for the lunar perigee it is [TT| [201 [2T| [221 [23] 

5w ~ 0.1 mas yr~ , (30) 



eq. (p9|) yields just 

s ~ 10"^ (31) 

in substantial agreement with the constraints of eq. (HD-eq. dSD obtained by Battat et al. [T2] 
as the outcome of a fit of modified models to LLR data. 

Finally, it may be interesting to look at the Earth and the LAGEOS satellites. Indeed, by 
looking at, say, the nodes of both LAGEOS and LAGEOS II, and the perigee of LAGEOS II 
it is possible to constrain s as previously done with the perihelia of three planets. Generally 
speaking, a major source of systematic uncertainty in the knowledge of the orbit of a terrestrial 
spacecraft is represented by the mismodeling in the first even zonal coefficient C2.0 accounting 
for the terrestrial quadrupole mass moment because of the resulting relatively huge secular 
precessions of the node and the perigee [2l]. By assuming [25] SC 2,0 ~ 10~^°, corresponding to 
uncertainties in the orbital elements of LAGEOS and LAGEOS II of the order of ~ 100 — 200 
mas yr~-^, it turns out 

Sx = -s°^ < 1.4 X 10-^ (32) 



s 



y 



-s°2 < 2 X 10-^ (33) 

s^ = -s03 < 1.3 X 10"^ (34) 

It can be noticed that eq. (I32p -eq. ()34p are neither competitive with the planetary bounds of 
eq. (f25|) -eq. pT)) nor with the lunar ones of eq. ([l])-eq. ([5]). 
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